Study on dynamic response and long-term settlement of silty soil around Shanghai metro tunnel

At present, the method for calculating long-term tunnel settlement predictions under metro loading considers only one working condition of passenger loading, which is inconsistent with actual working conditions. To establish a tunnel settlement model that accounts for variations in passenger flow, this study uses data mining methods to categorize metro operation into three working conditions: "peak period, secondary period, and low period." The impact of these passenger flow conditions on the dynamic response of the soil around the tunnel is analyzed. Then, based on the principles of calculus, a calculus-based prediction model is established to consider the changing patterns of metro passenger flow. The model is applied to analyze the long-term settlement characteristics of Shanghai Metro Line 10. The results indicate that, under identical conditions, soil displacement and dynamic deviatoric stress around the tunnel increase with passenger capacity. The calculus prediction model aligns more closely with actual working conditions than the conventional model. The predicted tunnel settlement of Shanghai Metro Line 10 after 20 years of operation is approximately 37.07 mm, with most settlement occurring in the early stages, primarily due to cumulative plastic deformation of the soil.

method 23 , dividing metro train operation into three working conditions: "peak period, secondary period, and low period."A three-dimensional coupling model of the metro-tunnel-soil system is established using Midas GTS NX software to analyze the dynamic response of the metro under each working condition.Then, based on the principles of calculus, an explicit calculus prediction model that considers the three operating conditions is established.The model's strengths and weaknesses are evaluated using Shanghai Metro Line 1 as an example.Finally, the model is applied to analyze the long-term settlement characteristics of the shield tunnel in the silty soil area of Shanghai Metro Line 10.

Model parameters
The tunnel for Shanghai Metro Line 10 has an outer diameter of 6.2 m, an inner diameter of 5.7 m, and is located at a depth of 12.5 m.It traverses soil layers ② 3a and ② 3b , with the physical parameters of these layers presented in Table 1.The tunnel segment is constructed of C55 concrete, while the track bed and floating slab are made of C30 concrete.The soil is simulated using the Mohr-Coulomb model, and the lining, roadbed, and floating slab are modeled using an elastic model.Welded contacts are employed between all model components.
Taking into consideration the soil environment surrounding the tunnel and the range of vibration influence 22 , the transverse width of the model (X) is set at 60 m (equivalent to 5 times the tunnel diameter from the tunnel's center to the outside), the vertical depth (Z) is 50 m, and the longitudinal length (Y) is 80 m.The accuracy of the dynamic analysis model's results is high when the mesh size is smaller than 1/14 of the shear wavelength (2.7 m) 24 .Thus, the grid size is set between 0.7 and 2 m to balance modelling efficiency, computational accuracy, and computational efficiency.For geostress equilibrium, a statically constrained boundary is used to zero out the displacements.To prevent wave reflection during the application of metro loads, an artificial viscoelastic boundary is utilized.The specific 3D numerical model is depicted in Fig. 1.

Analysis of metro passenger flow rules
The curve illustrating the change inmetro passenger flow is shown in Fig. 2 25 .It is evident from the figure that passenger flow fluctuates significantly within each time period.If only one operating condition were considered in the settlement calculation, the result would not align with the actual operating conditions.Consequently, this study utilizes data mining and the 3-4-5 rule 23 to segment metro operations based on passenger volume characteristics across various time periods.
Upon comparing the number of passengers for each condition, a significant discrepancy is observed.The number of passengers during the peak period is 1.93 times that of the secondary period and 3.92 times that of the low period.Hence, the variations in metro load due to passenger flow changes should be given due consideration.These fluctuations in passenger volume must also be factored into the related settlement calculations.

Simulation of metro cyclic load
Shanghai Metro Line 10 utilizes Type A Metro with a maximum axle weight of 16 tons, capable of accommodating 2,500 passengers at full capacity and traveling at a speed of 80 km/h 6 .Midas GTS NX software features a metro moving load function that automatically generates the moving load data based on the metro's parameters.The load is set in the model with the metro running from y = 0 m to y = 80 m, treating full capacity as the peak period load condition (with secondary and low periods scaled accordingly).The analysis focuses on the model's central section (y = 40 m), where the metro's dynamic load is depicted in Fig. 3, and the applied load is illustrated in Fig. 4.

Analysis of model results
The numerical results are analyzed under three metro loading conditions to investigate the variation rule of foundation settlement concerning the vertical and horizontal distance from the tunnel after the first loading, as well as the influence of operating conditions on foundation settlement and dynamic deviatoric stress.

Tunnel displacement characteristics after the first loading
Long-term tunnel settlement is primarily related to the soil deformation at the tunnel's base.Taking the soil at the bottom of the tunnel as the observation point, the time curve of the tunnel's vertical displacement for each working condition is depicted in Fig. 5.It is observed from Fig. 5 that the displacement of the soil beneath the tunnel for different working conditions follows a similar period, increasing and then decreasing with the metro's passage.The peak vertical displacement of the soil occurs when the fourth car of the metro runs through the observation section.The foundation deformation caused by metro vibrations increases with passenger volume.The maximum displacement during the peak, secondary, and low periods is 5.84 mm, 4.71 mm, and 4.05 mm, respectively.As passenger flow decreases from peak to secondary and then to low, the displacement decreases by 19.35% and 30.65%, respectively.It can be concluded that the effect of passenger volume changes cannot be neglected in the study of tunnel settlement.
The vertical displacements of the three working conditions were plotted as Fig. 6, which reveals the characteristics of displacement development in space.The figure shows that the vertical displacements for the three conditions exhibit the same spatial development trend.Displacement is greatest in the tunnel section and then gradually decreases horizontally outward and vertically downward from the tunnel as the center.This is because

Dynamic deviating stress response rule
The dynamic stress response level of the foundation soil is a critical factor affecting the long-term settlement and deformation of tunnels.The metro-induced dynamic deviatoric stress is typically represented by the second deviatoric stress invariant q d , which can be calculated from the soil stress components using Eq. 1 26 .
where: q d is the dynamic deviatoric stress; J 2 is the second stress invariant; σ xx , σ yy , σ zz , τ xy , τ xz , and τ yz are the stress components under the metro load.Figure 7 illustrates the dynamic deviatoric stress response curve of the soil beneath the tunnel.It can be seen that the dynamic deviatoric stress initially increases with depth, reaching a peak at 1.8 m, and then gradually decreases.This phenomenon is attributed to the superposition effect of vibrations caused by the different attenuation speeds of the body wave and the Rayleigh wave in the soil layer 27 , typically occurring 1-2 m below the tunnel 18,19 .The magnitude and depth of this effect are influenced by factors such as the type of soil layer, surrounding buildings, and tunnel structure.During actual construction, reinforcement measures should be implemented in the vibration amplification area to mitigate the impact of vibration amplification on settlement.
The peak value of the dynamic deviatoric stress decreases as passenger load decreases, with a 19.4% reduction from the peak period to the secondary period, followed by a 13.8% reduction from the secondary period to the low period.Therefore, when calculating the long-term settlement of the tunnel, the influence of passenger loading conditions must be considered to accurately account for the variation in dynamic response due to different operating conditions. (1)

Calculation model of long-term settlement of tunnel
To account for cumulative plastic deformation under cyclic loading and consolidation settlement due to the dissipation of cumulative pore water pressure, the hierarchical summation method is utilized to calculate the long-term settlement induced by metro vibration.

Calculation model of undrained cumulative deformation
The cumulative deformation law of a silty soil foundation under metro loading can be expressed by the exponential model in Eq. 2 28 .
where: ɛ p is the cumulative plastic deformation; q s is the static deviatoric stress; q d is the dynamic deviatoric stress amplitude; q f is the static destructive deviatoric stress; N is the number of cyclic load; a, m, t, and b are constants that can be determined through dynamic triaxial testing.

Calculation model of undrained excess pore pressure
The relationship between excess pore water pressure and the number of load cycles in silty soil under metro loading can be expressed as Eq. 3 29 .
where: u is the excess pore water pressure; q d is the dynamic deviatoric stress; q s is the static deviatoric stress; q f is the static destructive deviatoric stress; p c is the initial effective consolidation pressure of the soil body; α, β, γ, χ are the model parameters, which can be derived from the results of dynamic triaxial tests.

Calculation method of long-term cumulative settlement
The layered summation method is used for calculating tunnel settlement.The specific process is shown in Eqs.4-Eq.6 7 .
The settlement S d caused by the cumulative plastic deformation of the soil under metro loading can be calculated by the Eq. 4.
where: ɛ pr the cumulative plastic deformation of the soil in the rth layer; H r is the thickness of the soil in the rth layer; z is the number of soil layers in the compression layer.
The consolidation settlement S v due to dissipation of cumulative pore water pressure, which can be calculated by the cumulative pore pressure calculation model combined with the Terzaghi one-dimensional consolidation theory by the Eq. 5.
(2) where: u r the excess pore water pressure of the rth layer; mv r is the volume compression coefficient of the rth layer of soil; U r is the degree of consolidation of the rth layer of soil, taking 100%.The total sedimentation S Q is as the Eq. 6.

calculus prediction model
In the long-term operation of the metro, the variation in dynamic response between different working conditions cannot be ignored.To comprehensively consider the influence of each working condition, the conventional model of long-term settlement of underground tunnels is improved by incorporating the principles of calculus.Thus, a calculus-based model that accounts for changes in passenger flow is established.Assume the number of metro runs per day is A, with i runs during the peak period, j runs during the secondary period, and k runs during the low period.Integrate the number of load cycles N in Eq. 2 to divide the settlement caused by the first day of operation into three parts.As shown in Eq. 7.
where: h times is the number of loading times per vehicle, h = 12 21 ; ɛ p1I , ɛ p1J , ɛ p1K are the cumulative plastic deformations of peak, secondary, and peak period conditions in a day, respectively, which can be calculated by the Eqs.8-10.
The total deformation after nth day is ɛ pnQ , which can be expressed by Eq. 11.
where: ɛ pnIQ , ɛ pnJQ , ɛ pnKQ are the total amount of cumulative plastic deformation for peak, secondary, and peak period conditions after n days, respectively, which can be derived from the Eq. 12 ~ Eq. 14.
Likewise, the cumulative pore pressure after n days, considering the three operating conditions, can be calculated by Eq. 15.
In summary, the improved calculus prediction model can be formulated as Eq.16.
(5) where: S nQ , S nd , and S nv are the total settlement, cumulative plastic deformation settlement, and consolidation settlement with dissipation of excess pore pressure after day n, respectively.ɛ pnrQ is the cumulative plastic deformation of the rth layer of soil after day n; u nrQ is the excess pore water pressure in the rth layer after day n.

Reliability verification of calculus model
Figure 8 displays the predicted results of the tunnel settlement for Shanghai Metro Line 1, utilizing both the conventional model 30 and the improved model.The Mean Square Error (MSE) of the calculus model is observed to be 6.82, which is lower than the of the conventional model at 15.82.Consequently, the improved calculus model provides a more accurate prediction of the long-term settlement of metro tunnels compared to the conventional model.

Long-term settlement analysis of Shanghai Metro Line 10 tunnel
Dynamic triaxial test results from the silty soil surrounding the Shanghai Metro Line 10 tunnel, as reported in the literature [31], yield the following parameter values for Eqs. ( 2) and ( 3): a = 0.11, m = 1.42, t = 1.07, b = 0.19, c = 0.79, g = 1.45, k = 0.94, and l = 0.17.The first metro on Shanghai Metro Line 10 begins operation at 5:25 AM, and the last metro is at 10:30 PM, with metros running at approximately 5-min intervals.Utilizing the improved calculus prediction model, the long-term settlement 20a after the railway's opening has been forecasted, with the results illustrated in Figs. 9 and 10.As depicted in Fig. 9, the settlements of the Shanghai Metro Line 10 tunnel after 1, 5, 10, and 20a of operation are 24.49mm, 30.53 mm, 33.63 mm, and 37.07 mm, respectively.It is evident that the long-term settlement of tunnels in a silty soil layer, induced by metro train operations, increases over time.However, the majority of the settlement occurs within the first five years following the commencement of metro services, contributing to 82% of the total settlement.After this period, the rate of settlement decelerates and tends to stabilize.The settlement attributable to the cumulative plastic deformation of the soil amounts to 28.57 mm, representing 77% of the total settlement.Consequently, in relevant projects, it is crucial to prioritize disaster prevention and control related to tunnel settlement during the initial years of metro operation, with a particular focus on the cumulative plastic deformation of the soil.
After two decades of subterranean operation, the settlement of the soil at different depths beneath the tunnel is exhibited in Fig. 10.The maximum settlement displacement, occurring at a depth of 1.81 m (the peak point of dynamic deviatoric stress) beneath the tunnel, is 13.14 mm.This indicates that the long-term settlement of the soil in the 0 to 1.81 m range constitutes 23.93 mm, accounting for 64.55% of the total settlement.Therefore, the deformation of soil due to tunnel settlement predominantly takes place within the 1.81 m depth below the tunnel, specifically in the area where dynamic deviatoric stress vibration is amplified.In relevant projects, reinforcement measures should be concentrated at this depth to mitigate the detrimental effects of differential settlement on the tunnel structure and its operation.

Conclusions
In this study, the effect of changes in metro passenger flow on the long-term settlement of tunnels was considered.The operation of the underground is divided into three working conditions.Based on this classification, a three-dimensional numerical model of the metro-tunnel-strata interaction is established to compare and analyze the dynamic stress response under these conditions.Subsequently, drawing on calculus principles, an explicit calculus prediction model is formulated, which incorporates changes in passenger flow and investigates the settlement development law of tunnels in silty soil.The main conclusions are as follows: 1.The data mining method is utilized to categorize metro operations into three working conditions: peak, secondary, and low peak periods.Significant differences in passenger load among these conditions are identified.In related projects, it is important to account for the variations in metro load due to passenger flow fluctuations and to consider the impact of these changes in related research.2. The soil displacement around the tunnel due to metro operations decreases with increasing distance from the tunnel.The dynamic deviatoric stress shows an amplification zone below the tunnel, where the stress increases with depth, peaks at 1.81 m, and then decreases gradually.The vertical displacement and dynamic deviatoric stress response of the soil surrounding the tunnel are exacerbated as passenger capacity increases.3. The explicit long-term settlement model for metro tunnels is refined using calculus principles, and a calculus model that accounts for the three working conditions is developed.The reliability of this calculus model is confirmed using the measured data from Shanghai Metro Line 1.The results indicate that the improved calculus prediction model is more appropriate for studying long-term tunnel settlement.4. The total tunnel settlement for Shanghai Metro Line 10 after 20 years is 37.07 mm, with 82% occurring within the first 5 years of the metro's operation.The settlement resulting from the cumulative plastic deformation of the foundation soil constitutes 77% of the total settlement, which is the primary component of long-term tunnel settlement.Tunnel settlement is predominantly caused by soil deformation in the vibration amplification zone, with settlement from this zone accounting for 64.55% of the total.In related projects, it www.nature.com/scientificreports/ is essential to implement reinforcement measures within the vibration amplification zone.Attention should be directed toward disaster prevention and control of tunnel settlement during the metro's preoperational period, with a particular focus on the cumulative plastic deformation of the soil to prevent uneven settlement from adversely affecting the tunnel's structure and operation.
This study builds on existing research on the prediction method for long-term tunnel settlement by incorporating the impact of passenger flow changes, establishing a settlement calculation model that considers these changes.However, it is important to note that the calculation only takes into account the static load from the structure itself and the effects of metro loads, and that tunnel settlement due to metro operation represents only a part of the total long-term settlement.In the long-term operation of the metro, factors such as construction activity around the tunnel, soil rheology, and changes in the water table will also impact tunnel settlement, possibly leading to deviations between predicted and actual settlements.Therefore, a comprehensive calculation method for tunnel settlement that considers various factors still requires in-depth study.

Figure 2 .
Figure 2. Metro passenger flow curve over time.

Figure 7 .
Figure 7. Dynamic deviatoric stress response of the soil underneath the tunnel.

Table
Table of modelled soil parameters.